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Abstract 

Multifractal  characteristics  of  the  upper  layer  thermal  structure  in  the  southwestern  Greenland  Sea,  Iceland  Sea,  and 
Norwegian  Sea  (GIN  Sea)  are  analyzed  using  high-resolution,  digital  thermistor  chain  data.  The  energy  spectrum  at  20 
m  depth  (cold  sublayer)  shows  the  existence  of  a  spike  at  the  scale  of  approximate  3  km  representing  the  chimney  scale. 
The  graph  dimension  varies  from  higher  values  such  as  1.89  at  the  surface  to  1.44-1.50  in  the  warm  intermediate  layer. 
The  stationarity  decreases  from  the  ocean  surface  to  the  warm  intermediate  layer.  However,  the  information  dimension 
varies  slightly  (0.92  to  0.90)  that  indicates  low  singularity. 
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1.  Introduction 

The  Greenland  Sea,  Iceland  Sea,  and  Norwegian  Sea  (GIN  Sea)  are  key  regions  in  the  advective-convective  system 
with  various  stages  of  modification  that  links  the  polar  ocean  with  the  North  Atlantic  (Fig.  1).  Because  of  the  import 
and  modification  of  water  masses  a  large  number  of  regional  water  types  can  be  encountered.  The  North  Atlantic  Water 
(NAW)  is  relatively  wanner  and  saline  (T  >  2  °C,  S  >  34.9  ppt).  The  Arctic  Water  (AW)  is  cooler  and  fresher  (T  <  0 
°C,  S  <  34.7  ppt)  [1],  Different  water  masses  encountered  in  the  GIN  Sea  interface  and  form  fronts  and  eddies  that  not 
only  determine  water  bodies  with  different  hydrographic  characteristics  but  also  the  regional  biological  systems.  The 
major  circulation  pattern  of  the  GIN  Sea  is  characterized  by  the  Northward  flowing  Norwegian  Atlantic  Current 
(NAC)  in  the  east  and  the  southwestward  flowing  East  Greenland  Current  (EGC)  along  the  East  Greenland  shelf  and 
leaving  the  Iceland  Sea  via  the  Denmark  Strait. 

Convection  in  GIN  Sea  is  seen  as  a  globally  important  process  in  which  air-sea  interactions  influence  oceanic 
circulation  through  the  production  and  ventilation  of  deep  and  intermediate  waters.  The  key  dynamic  elements  of 
oceanic  convection  are  taken  to  be  the  individual  plumes,  clusters  of  plumes  or  called  chimneys,  and  eddies  that  are  the 
consequence  of  chimneys  aging  in  a  rotating  frame.  Two  major  features,  nonstationarity  and  intermittency,  should  be 
first  investigated  in  order  to  understand  oceanic  convective  process  or  the  secondary  circulation  across  oceanic  fronts 
[2].  Question  arises:  How  can  we  determine  upper  ocean  nonstationarity  and  intermittency  from  observational  data? 
This  paper  describes  a  multi-fractal  analysis  on  a  high-resolution  temperature  dataset  to  obtain  the  nonstationarity  and 
intermittency  of  the  upper  layer  (300  m  depth)  in  the  southwestern  GIN  Sea. 


2.  Thermistor  chain  data 

In  July-August  1987,  fine-resolution  temperature  data  are  collected  on  board  of  M /V  SEA  SEACHER  by  the  Royal 
Navy’s  Admiralty  Research  Establishment  (ARE)  using  a  digital  thermistor  chain  (280  m  long,  100  sensor  pods)  with  a 
single  200  km  straight-line  tow  (4  knots  speed)  near  69  °N,  18  °W  (Fig.  2)  in  the  east  edge  of  EGC  [3].  The  upper  ocean 
(surface  to  280  m  depth)  was  sampled  every  0.9  s,  obtaining  a  temperature  profile  about  every  2  m  along  the  track.  Each 
sensor  pod  of  the  chain  measures  temperature,  and  about  one  in  five  also  measures  pressure,  allowing  the  depth 
distribution  of  temperature  to  be  deduced. 
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Fig.  1 .  Schematic  representation  of  bathymetry  and  currents  in  the  GIN  Sea. 


Fig.  2.  Track  along  which  the  thermistor  chain  data  were  taken.  The  dotted  curve  represents  1000  m  depth  contour  (after  Scott  and 
Killworth  [3]). 


Fig.  3  shows  contour  plot  of  temperature  ranging  from  2  to  8  °C  with  0.2  °C  increment  on  a  vertical  cross-section 
between  the  two  marked  locations  ‘ b ’  and  ‘ e ’  (Fig.  2).  As  pointed  by  Scott  and  Killworth  [3],  the  temperature  shows  a 
small-scale  variability  with  highly  irregular  nature.  The  total  length  of  the  temperature  cross-section  is  90  km  with  each 
interval  of  6  km.  The  thin  surface  layer  (depth  around  5  m)  temperature  is  about  +7.5  °C.  Below  the  surface  layer,  there 
is  a  first  thermocline  with  temperature  rapidly  decreasing  with  depth  to  25  m.  The  vertical  gradient  in  the  first  ther- 
mocline  is  very  strong  (~0.45  °C/m).  Cold,  relatively  uniform  sublayer  exists  below  the  first  thermocline  from  25  to  70 
m,  reaching  a  minimum  temperature  of  -1.5  °C.  Below  the  cold  sublayer,  there  is  a  second  thermocline  (70-90  m, 
thickness  around  20  m)  where  the  temperature  increases  with  depth  with  a  vertical  temperature  gradient  around  0.135 
°C/m.  Below  the  second  thermocline,  there  exists  a  warm  intermediate  layer  (~  100-300  m)  with  a  maximum  temper¬ 
ature  of  +1.2  °C.  Note  that  the  thermal  characteristics  are  different  between  the  two  thermoclines.  The  first  (second) 
thermocline  features  the  decrease  (increases)  of  temperature  with  depth.  The  chimneys  are  observed  below  the  char¬ 
acteristic  cold  sublayer  (70  m  depth)  with  width  of  3  km  (Fig.  4).  The  water  is  about  1.2  °C  cooler  within  than  out  of  the 
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Fig.  3.  Temperature  cross-section  obtained  from  the  thermistor  chain  data  collected  along  the  track  from  Station-b  to  Station-e  (after 
Scott  and  Killworth  [3]). 
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Fig.  4.  An  expansion  of  central  part  of  Fig.  3  showing  30  km  section  between  Station-c  and  Station-d  (after  Scott  and  Killworth  [3]). 


chimneys.  The  chimneys  reach  the  maximum  thermistor  chain  depth  (~280  m).  Containing  water  of  temperature  down 
to  about  0  °C,  the  chimneys  appear  to  punch  the  cold  water  clearly  through  the  warm  intermediate  layer  water  with  a 
maximum  temperature  of  +1.2  °C. 


3.  Power  spectra 

What  is  the  inherent  thermal  variability  identified  from  this  high-resolution  temperature  data  with  multi-layer  and 
chimney  structures?  What  are  the  statistical  properties?  Before  answering  these  questions,  we  should  first  investigate 
the  stationarity  of  the  temperature  field.  For  a  given  depth,  the  temperature  data  is  a  function  of  the  horizontal 
coordinates,  x, 

T,  =  T(Xi),  x,  =  il,  i  =  0,l,...,A,  A  =  L/l,  (1) 

where  /  =  2  m,  is  the  horizontal  resolution  of  the  data,  and  L  is  the  total  horizontal  scale  of  the  data  set.  Our  data  set 
has  24,000  temperature  profiles,  that  is,  A  -  24,000. 

Spectral  analyses  of  temperature  field, 

Ej  =  E(kj),  kj=j/L ,  y  =  1,2,...,  A/2.  (2) 

at  all  depths  were  done,  but  for  the  sake  of  brevity  and  to  elucidate  the  important  points,  only  spectrum  at  20  m  depth 
(cold  sublayer,  Fig.  5)  is  shown.  A  Bartlett  window  was  used  to  taper  the  ends  of  each  series  before  calculating  the 
power  spectra  to  reduce  the  spectral  leakage  in  the  wavenumber  domain.  Fig.  5  shows  the  existence  of  a  spike  on 
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Fig.  5.  Power  spectrum  of  temperature  filed  at  20  m  depth. 


log2(A£/2)  ~  6.5, 

corresponding  to  a  scale  of  approximate  3  km  in  the  thermal  variability.  This  length  scale  (3  km)  coincides  with  the 
chimney  scale,  which  may  imply  the  existence  of  a  linkage  between  the  cold  sublayer  and  the  chimney  formation. 


4.  Stationarity 

For  a  scaling  process,  one  expects  power  law  behavior  [4], 

E{k)  oc  k~f,  (3) 

over  a  large  range  of  wavenumber  k.  The  spectral  exponent  /S  contains  information  about  the  degree  of  stationarity  of 
the  data  [5,6],  If  0  <  1,  the  field  is  stationary;  if  1  <  f)  <  3,  the  field  contains  nonstationary  signal  with  stationary  in¬ 
crements  and  in  particular,  the  small-scale  gradient  field  is  stationary;  if  /?  >  3,  the  field  is  nonstationary  with  non¬ 
stationary  increments. 

The  power  spectra  for  all  the  depths  have  multi-scale  characteristics  with  the  spectral  exponent  /?  in  the  range  of 

1  <  P  <  2,  (4) 

which  means  the  temperature  field  of  the  southwestern  GIN  Sea  sublayer  is  nonstationary  with  stationary  increments 
(Fig.  5). 


5.  Structure  functions 

Since  the  thermistor  chain  data  set  has  stationary  increments,  we  should  study  the  statistical  characteristics  of  the 
gradient  field, 

|AT(r;x)|  =  |T(xi+r)  -  F(jc,)|,  /  =  0, 1, . . . ,  A  -  r, 

where  r  denotes  the  lag  between  two  data  points.  Obviously,  r  is  inversely  proportional  to  the  wavenumber  k, 

1 


(5) 
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Structure  functions  are  used  to  depict  the  thermal  variability.  The  gth-order  structure  function  is  defined  by  the 
mean  of  the  gth  power  of  the  gradient  field  |A7’(r;x)|, 

S(r,q)  =  <|Ar(r;*)H  I Wcx)?.  (6) 

/A  .  /« 0 

For  example,  r  —  1,  q  =  1,  the  structure  function 

=  -T3T  1^/+')  ~  r(x,)|,  (7) 

71  1  0 

represents  the  average  magnitude  of  gradient. 

Near-linear  dependence  of  log2[S(r,g)]  on  log2(r)  is  found  with  different  ^-values  from  0.5  to  4.0  (Fig.  6).  The 
straight  lines  with  different  slopes  show  that  the  structure  functions  for  the  upper  ocean  temperature  in  the  southwestern 
GIN  Sea  satisfies  the  power  law 

S(r,q)  oc^(,),  (8) 

with  the  exponent  £(q)  depending  on  q.  Since  S(r,  0)  =  1,  the  power  £(g)  should  be  0  for  q  =  0.  Our  computation  (Fig.  6) 
agrees  quite  well  with  earlier  studies  [7,8]:  Z(q)  is  monotonically  and  near-linearly  increasing  with  q.  Thus,  it  may  be 
represented  by 

Z(q)=H(q)q,  (9) 

where  H(q)  is  nearly  a  constant. 

The  power  of  the  structure  function,  Z(q),  is  computed  for  the  whole  dataset.  Fig.  7  shows  the  dependence  of  £(<?)  on 
q  for  selected  levels.  All  the  curves  (near-linear)  converges  at  q  =  0  [£(0)  =  0],  and  show  three  different  patterns:  (a)  slow 
increasing  with  q  (surface  layer,  0  m),  (b)  intermediate-rate  increasing  with  q  (second  thermocline,  80  m),  (c)  fast  in¬ 
creasing  with  q  (first  thermocline,  20  m;  cold  sublayer,  40-60  m;  intermediate  warm  layer,  below  100  m). 

The  structure  function  for  q  =  1,  S(r,  1),  is  often  used  to  determine  the  statistical  characteristics  of  the  data  such  as 
stochastically  continuous  and  stationarity.  When  the  structure  function  S(r,  1)  does  not  depend  on  r 

S(r,  1)  =  const,  (10) 

which  denotes  exact  stationary.  Eq.  (10)  is  equivalent  to 

f(l)=#(l)=0. 
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Fig.  7.  Dependence  of  the  structure  function’s  power,  £(■?),  on  q  and  depth. 


Table  1 

Power  of  the  structure  function  S(r ,  1)  and  the  dimension  of  the  GIN  Sea  thermal  field  T 


Depth  (m) 

0 

20 

40 

60 

80 

100 

120 

140 

200 

C(l)  =  *. 

0.11 

0.43 

0.47 

MEM 

0.40 

0.56 

Ds(T) 

1.89 

1.57 

1.53 

KS 

1.60 

1.44 

mem 

Thus,  H\  =H(l),  represents  statistical  characteristics  of  data.  Since  the  power  ((1)  (i.e.,  //(I))  varies  with  depth 
(Table  1)  between  0.11  at  the  surface  and  0.56  at  120  m  depth  (intermediate  warm  layer),  the  thermal  field  in  the 
southwestern  GIN  Sea  is  nearly  stationary. 

If  g[T)  is  defined  as  the  graph  of  T(x),  whose  dimension  can  be  defined  by  [9] 

Dg(T)  —2  —  H\.  (11) 

For  a  stochastically  continuous  (H\  =  1)  data  T(x),  the  graph  should  be  a  smooth  curve,  whose  dimension,  Dg(T) 
should  be  1.  If  the  graphs  g(T)  fill  the  whole  space  (exact  stationary), 

Dun  =  2,  (12) 

which  corresponds  to  H\  =0. 

Table  1  shows  that  the  dimension  of  the  southwestern  GIN  Sea  upper  thermal  field  varies  (multi-dimension 
structure)  from  higher  values  such  as  1.89  (surface:  0  m),  1.78  (second  thermocline:  80  m)  to  lower  values  such  as  1.57 
(first  thermocline:  20  m),  1.52-1.53  (cold  sublayer:  40,  60  m),  and  1.44-1.50  (warm  intermediate  layer:  120, 140, 200  m). 
The  decreasing  order  of  stationarity  is:  ocean  surface,  second  thermocline,  first  thermocline,  cold  sublayer,  and  warm 
intermediate  layer. 


6.  Singular  measures 

Normalized  small  scale  absolute  gradient 


(13) 
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is  used  to  identify  the  intermittency  of  the  thermal  field.  The  running  average  of  r  normalized  values  are  computed  by 

e{r;x,)=-  *0;*/).  i  =  0, 1, . ..  ,A  -  r. 


The  mean  of  the  qXh  power  of  s(r;;c,) 

M{r,q)  =  (s(r,Xj)9)  =  -±-  £>(r; 

/l  r  (=0 

is  defined  as  the  ^th-order  singular  measure.  Obviously,  for  q  =  0, 

M{r,  0)  =  1. 

For  q  ~ 

M(r,  1)  =  (e(r;x,))  =  -i-  ^[e(r;x,)]  =  — l—  £  ;  £  e(l;*/) 


/'=*0 


(14) 


(15) 


(16) 


(17) 


The  singular  measures  are  computed  for  all  depths.  For  simplicity,  M(r,  q)  for  the  depth  of  40  m  is  given  here  (Fig. 
8).  Near-linear  dependence  of  log2[Af (r,  g)]  on  log2(r)  is  found  with  different  q- values  from  0.5  to  4.0.  The  straight  lines 
with  different  slopes  show  that  the  singular  measures  with  various  q  for  the  upper  layer  temperature  in  the  southwestern 
GIN  Sea  satisfies  the  power  law 


M(r,q)<xr-KM,  q>  0, 

with  the  power  K(q)  varying  with  q.  From  Eqs.  (16)  and  (17),  we  have 
AT(0)  =^(1)  =  0. 

Several  characteristics  are  found  from  Fig.  9:  The  power  K(q)  is  a  convex  function 

^>0 

dq2  ’ 

for  all  q  and 

K(q)  <  0  only  if  0  <  q  <  1 


(18) 

(19) 

(20) 

(21) 


Depth  40m 


Fig.  8.  Singular  measures  for  different  ^-values  (depth  =  40  m). 
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q 


Fig.  9.  Dependence  of  the  singular  measure’s  power,  K(q),  on  q  and  depth. 


which  reflects  the  fact  that,  in  this  range,  taking  a  gth  power  necessarily  reduces  the  fluctuation  of  e(r,x,);  and  otherwise 


K(q)  ^  0,  if  q  >  0. 


(22) 


Following  [10,11]  we  may  define  a  function 


C(q)  = 


K(q) 
q- 1' 


(23) 


For  q  — >  1,  we  use  FHosptaT s  rule  to  define  a  straightforward  measure  of  inhomogeneity  in  the  sense  of  singular 
measure  [4]: 

C,=C(1)=*'(1)>0,  (24) 

which  is  called  the  intermittency  parameter.  The  larger  the  value  of  C\  the  lager  the  intermittency  and  singularity  the 
data  set  has.  The  intermittency  parameter  C \  varies  from  0.08  to  0.1. 


7.  Mean  multi-fractal  plane 

Statistical  moments  of  all  orders  (gth  power)  can  be  computed,  however,  the  first  moment  (q  =  1)  provides  im¬ 
portant  information  about  the  data.  The  parameter  C|  measures  the  degree  of  intermittency  in  the  system,  while  H\ 
measures  its  degree  of  nonstationarity.  The  plot  of  C\  versus  H\  called  the  mean  multi-fractal  plane,  shows  the  degree  of 
nonstaionarity  and  intermittency  (Fig.  10). 

Both  parameters  have  geometrical  interpretations  as  co-dimensions:  information  dimension  and  graph  dimension 
[4],  The  information  dimension,  represented  by  (1  —  Ci),  is  a  first-order  estimate  of  sparseness  of  strong  gradient 
distributed  in  the  system.  The  graph  dimension,  represented  by  (2  -  Hi),  is  a  first-order  estimate  of  roughness  in  the 
system.  Both  parameters  have  analytical  meanings:  Ct  (sparseness)  is  related  to  singularity  and  H\  (roughness)  is 
connected  to  the  lack  of  differentiability.  The  southwestern  GIN  Sea  upper  layer  thermal  field  shows  multi-fractal 
characteristics  (Fig.  10)  in  terms  of  Ci  and  H\  We  find  larger  variation  in  H\  (0.11-0.56)  than  in  C\  (0.08-0.10). 
Identification  of  {H\,C\)  ranges  helps  to  select  realistic  model  to  describe  the  field  [4]. 
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Fig.  10.  Mean  multi-fractal  plane  for  the  southwestern  GIN  Sea  thermal  field. 


8.  Conclusions 

The  multi-fractal  analysis  provides  a  useful  framework  for  analyzing  ocean  data  when  complex  nonlinear  processes. 
The  upper  layer  thermal  structure  in  the  southwestern  GIN  Sea  has  the  following  features: 

(1)  The  energy  spectrum  at  20  m  depth  (cold  sublayer)  shows  the  existence  of  a  spike  at  the  scale  of  approximate  3  km 
in  the  thermal  variability.  This  length  scale  (3  km)  coincides  with  the  chimney  scale,  which  may  imply  the  existence 
of  a  linkage  between  the  cold  sublayer  and  the  chimney  formation. 

(2)  The  power  spectra  for  all  the  depths  have  multi-scale  characteristics  with  the  spectral  exponent  P  in  the  range  of  (1, 
2).  This  means  the  temperature  field  of  the  southwestern  GIN  Sea  sublayer  is  nonstationary  with  stationary  incre¬ 
ments. 

(3)  The  structure  function  has  multi-fractal  characteristics,  i.e.,  the  power  of  the  gth-order  structure  function  is  mono- 
tonically  and  near-linearly  increasing  with  q.  However,  the  rate  of  such  an  increasing  varies  with  depth:  (a)  slow 
increasing  with  q  (surface  layer,  0  m),  (b)  intermediate-rate  increasing  with  q  (second  thermocline,  80  m),  (c)  fast 
increasing  with  q  (first  thermocline,  20  m;  cold  sublayer,  40-60  m;  intermediate  warm  layer,  below  100  m). 

(4)  Two  parameters  for  the  first  moment  (q  =  1)  structure  function  and  sigular  measure,  H\  and  C \  are  important  to 
represent  the  multi-fractal  characteristics  in  terms  of  co-dimensions  (information  and  graph).  The  graph  dimension, 
represented  by  (2  —  Hi),  is  a  first-order  estimate  of  roughness  in  the  system.  The  information  dimension,  represented 
by  (1  -  Ci),  is  a  first-order  estimate  of  sparseness  of  strong  gradient  distributed  in  the  system.  Both  parameters  have 
analytical  meanings:  Ci  (sparseness)  is  related  to  singularity  and  H\  (roughness)  is  connected  to  the  lack  of  differ¬ 
entiability. 

The  graph  dimension  of  the  southwestern  GIN  Sea  upper  thermal  field  varies  (multi-dimension  structure)  from 
higher  values  such  as  1.89  (surface:  0  m),  1.78  (second  thermocline:  80  m)  to  lower  values  such  as  1.57  (first  thermocline: 
20  m),  1.52-1.53  (cold  sublayer:  40,  60  m),  and  1.44-1.50  (warm  intermediate  layer:  120,  140,  200  m).  The  decreasing 
order  of  the  stationarity  is:  the  ocean  surface,  second  thermocline,  first  thermocline,  cold  sublayer,  and  warm  inter¬ 
mediate  layer.  However,  the  information  dimension  of  the  southwestern  GIN  Sea  upper  thermal  field  varies  slightly 
from  0.92  to  0.90.  This  indicates  that  singularity  is  very  low. 
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